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ABSTRACT 

0^ ' Context. In the course of the life of a massive star, wind- wind interaction can give rise to the formation of circumstellar nebulae which 

are both predicted and observed in the nature. 

Aims. We present generic model calculations to predict the properties of such nebulae for blue supergiants. 

Methods. From stellar evolution calculations including rotation, we obtain the time dependence of the stellar wind properties and of 

the stellar radiation field. These are used as input for hydro-calculations of the circumstellar medium throughout the star's life. 

Results. Here, we present the results for a rapidly rotating 12 Mq single star. This star undergoes a blue loop during its post main 
O ■ sequence evolution, at the onset of which its contraction spins it up close to critical rotation. Due to the consequent anisotropic mass 

$H ' loss, the blue supergiant wind sweeps up the preceding slow wind into an hour glass structure. Its collision with the previously formed 

0^ ' spherical red supergiant wind shell forms a short-lived luminous nebula consisting of two polar caps and a central inner ring. With 

j^ , time, the polar caps evolve into mid-latitude rings which gradually move toward the equatorial plane while the central ring is fading. 

These structures are reminiscent to the observed nebulae around the blue supergiant Sher 25 and the progenitor of SN 1987A. 

Conclusions. The simple model of an hour glass colliding with a spherical shell retrieves most of the intriguing nebula geometries 

discovered around blue supergiants, and suggests them to form an evolutionary sequence. Our results indicate that binarity is not 
^ ' required to obtain them. 

a^: 
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t~^ 1. Introduction may be difficult to observe, their violent interaction may release 

^^ enough energy to provide an observable nebula. 

OO During the course of their evolution, massive stars have strong Despite this simple and intriguing expectation, there are 

O winds which eject matter into their suiToundings. During their ^o far only few attempts to obtain quantitative prediction for 

L; ■ post-main sequence evolution, these stars can move back and the outcome of the described shell interaction ( see Bl ondin et 

.,-j forth from the blue to the red side of the Hertzsprung-Russell al. 1993, Martin et al. 1995 Podsiadlowski et aL 120051) . Within 

^ (HR) diagram and back to the red, with Uttle time spent an effort to describe this phenomenon through generic calcula- 

%-^ at intermediate effective temperatures (e.g., Langer 1991b). tions, which use detailed stellar evolution models as input for the 

Hydrodynamic considerations imply that each such transition circumstellar hydrodynamic modeling (Chija et. al., in prepara- 

does produce a circumstellar shell: When the star moves from tion), we focus here on the results for a rotating 12 Mo single 

the blue to the red side of the HR diagram, the slow red su- gjar 
pergiant (RSG) wind will be stalled by the high pressure of the 
previously created hot wind bubble, and will accumulate into a 

shell at the location where this p ressure equals the RSG wind 2. Computational method 
ram pressure (Garcia-Segura et al. 1 1996b] ). We call such a more 

or less stationary shell the RSG shell. When the star moves from As input for our circumstellar hydrodynamic calculations, we 

the red to the blue side of the HR diagram, the wind speed in- use the results of a stellar evolution calculation for a star of 

creases and the blue supergiant (BSG) wind plows up the pre- 12 Mq and a metallicity of Z = 0. 02. Sp ecifically, we utihze 

ceding RSG wind into a rapidly expanding shell, which we call Model F12B from Heger & Langer (I2000I I. which has an initial 

the BSG shell. rotational velocity of 328 km s ' . The code used to compute this 

Consequendy, we expect a spectacular circumstellar phe- "^o^el includes OPAL opacities, detailed nuclear networks, mass 

nomenon for stars undergoing so called blue loops, namely that l^^s according to Nieuwenhuijzen & Jager ( 1990), the physics of 

it triggers the formation of an expanding BSG shell, which rotation for the stellar interior, and rotationally modulated stellar 

will at some point smash into the previously formed stationary ^^^s, as described in Heger, Langer & Woosley ([2003. 

RSG shell. While both, the RSG and the BSG shell by itself. The evolution of the stellar model in the HR diagram is show 



a 



in Fig. lA.ll At core-H exhaustion, it moves to the RSG regime 
where it remains for 825 000 yrs (~60 % of the core-He burning 
Send offprint requests to: S. M. Chifa, s.m.chita@uu.nl fife time), before it undergoes a blue loop. It then stays in the 
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Table 1. Ejected mass (AM), momentum (Ap) and kinetic en- 
ergy (AE) during the various evolutionary phases of our stellar 
model. The evolutionary phase is identified in the first column: 
main sequence phase (MS), first red supergiant phase (RSGI), 
phase of rapid rotation (RR), blue supergiant stage (BSG), and 
second red supergiant phase (RSGII). The approximate duration 
of each phase is given in the second column. 



Phase 


At 


AM 


Ap 


AE 




lO'yr 


Mo 


lO^Sgcms-' 


10^5 erg 


MS 


19200 


0.43 


396 


1480 


RSGI 


825 


0.33 


71 


38 


RR 


25 


0.02 


7.2 


6.0 


BSG 


550 


0.11 


52 


68 


RSGII 


225 


0.13 


25 


12 



BSG regime of the HR diagram for the remaining part of core 
helium burning, before it moves back to the RSG regime where 
it explodes as a Type II supernova. 

As shown by Heger & Langer ( |1998l l, as the convective en- 
velope retreats during the onset of the blue loop, all its angu- 
lar momentum is concentrated in a small amount of mass in 
the top layers of the star by the time convection vanishes. Blue 
loops therefore provide a natural way to bring the stellar sur- 
face to close to critical rotation. This does also happen in our 
chosen stellar model (Fig. IA.2b . The limit of critical rotation is 
reached during the red-blue transition, which produces a brief 
period of strong, slow and anisotropic mass loss (Table 1). The 
strong mass loss then reduces the rotation rate of the stellar sur- 
face (Langer 19981 1, and the star settles at a rotation velocity of 
about 50kms ' in the BSG regime. 

To simulate the evolution of the circumstellar matter (CSM) 
we use the ZEUS 3D code developed by Stone & Norman ( |1992| l 
and Clark (1996 ). ZEUS 3D is an explicit non-conservative code 
that solves the hydrodynamical equations as partial, finite dif- 
ference equations on a fixed, staggered mesh grid. Radiatively 
optically-thin cooling is included by solving the energy equation 
implicitly according to Mac Low et al. ( I1989I I, and by using the 
plasma cooling curve of MacDonald & Bailey (11981) . We com- 
pute the evolution of the CSM during the main sequence and the 
early RSG stage in ID, with 4500 grid points over a radius of 
45 pc, and we assume an interstellar medium density of 1 cm"-'. 
After 100 000 yr into the first RSG stage, we map the ID model 
onto a 2D spherical grid to compute its further evolution. The 
inflow inner boundary condition is applied at 0.025 pc, and the 
outer boundary remained at 45 pc. The radial component of the 
grid is resolved with 1000 grid points, where 900 grid points 
are used for the inner 5 pc, and 100 grid points for the outer 
40 pc. The angular coordinate of 90 degrees is resolved with 
200 grid points. The method used here was applied before by 
Garcia-Segura et al. ( 119961 |1996a| and |1996bl l and van Marie et 
al. ( i2005 and 2007l l. 

We are using the time dependent mass loss rate and the ter- 
minal wind speed from the stellar evolution model as input in 
our central mesh point for the hydrodynamic calculations. The 
wind speed is obtained from the stellar escape velocity using the 
scaling proposed by Eldridge (120051). The wind anisotropy is de- 
scribed using the equations of Bjorkman & Cassinelli (|1993l l, 
as in Langer et al. (T999). For near-critically rotating stars, this 
provides a slow and dense equatorial outflow, and a fast wind 
in polar directions. We note that while the Bjorkman-Cassinelli 
mechanism has been criticized in the context of line driven winds 



(Owocki et al. 119961 ). it is unclear whether line driving does play 
a major role in the situation of near-critical rotation. 

The eff'ect of photoionization was included in the simulations 
by calculating the Stromgren radius along each radial grid line 
as described in Garcia-Segura et al. (119991 ) and van Marie et 
al. (120051 120071 and |2008 ). The number of ionizing photons 
is computed according to the effective temperatures and surface 
gravities of the stellar evolution model, by interpolating in a grid 
of model atmospheres for massive OB stars of solar metallicity 
computed with the FASTWI ND non -LTE code (Puis et al. |2005l ) 
as described in Lefever et al. (120071 ) . 

3. Results 

During its main sequence phase our 12 Mq star creates a hot bub- 
ble in the interstellar medium which, at core hydrogen exhaus- 
tion, is characterized by a radius of 30 pc and an internal energy 
density of 10 '^ergcm""*. Once the star has become a RSG, a 
slow (~ 50kms"'), dense and isotropic wind is injected into the 
computational domain (Fig. IA.2I ). This RSG wind accumulates 
at a distance of ~ 1.5 pc where its ram pressure is balanced by the 
hot bubble pressure, and forms a RSG shell (cf. Garcia-Segura et 
al. 1996b). At the end of the first RSG phase, this shell contains 
about 0.26 Mq. It is rather extended (Ar ^ 1 pc), and its central 
parts have condensed due to cooling. 

At the onset of the blue loop, the central star reaches close- 
to-critical rotation, and ejects a dense equatorial disk (Heger 
& Langer 1998). While mass and time scales differ, this phe- 
nomenon occurs quite analogous to the simulation of the out- 
burst of 77 Carinae by Langer et al. ( 119991 ). Like in this case, the 
ensuing BSG wind sweeps up the preceding slow wind material 
into an "hour glass" structure (Fig. |A.4t . On a time scale of a 
few lO'* yr, this hour glass expands into the sphere defined by the 
RSG shell, with a maximum velocity of ~ 130 kms"' (Fig. |A.3l ). 
The faster polar parts of the hour glass hit the inner edge of the 
RSG shell first. The collision creates a hot {T ^ 10^ K) and 
dense (n ^ lOcm"^) pair of polar caps. As time proceeds, the 
collision zone moves to lower latitudes of the RSG shell and be- 
comes more confined in latitude. At the same time, the interac- 
tion of the BSG wind with the equatorial disk defines a second, 
ring shaped collision zone in the equatorial plane, which expands 
with time with a velocity of 18 km s"' . 

Figure \T\ shows snapshots of the emissivity map, according 
to the employed cooling curve in our hydro simulations, for three 
slices in time, along with projection maps constructed from ro- 
tationally symmetric 3D-structures obtained from the 2D maps. 
Here, only emission from gas in the temperature range between 
lO"* K and 10^ K is considered, which is the dominant compo- 
nent. Hotter gas, which is formed from the reverse shock of the 
collision, might be observable in the X-ray regime; the peak lu- 
minosity of this component in our model is 10^^ erg s"' , which 
is achieved about 50 000 yr after the onset of the collision. At 
an early interaction stage, the radiation is dominated by two po- 
lar caps and one equatorial ring, later on by two mid-latitude 
rings and one fading smaller equatorial ring, and finally two 
mid-latitude rings at rather low latitude are visible. Those two 
rings gradually move to the equatorial plane while fading. The 
full time dependence of the emission structure is shown in an 
accompanying movie which is available from the A&A website. 

The energy budget for the collision of the polar caps of 
the hour glass with the RSG shell follows directly from the 
stellar properties. The polar caps have an emissivity of / ^ 
10"^' erg cm"-' s"' in a volume of V ^ 47rr^Ar = 4 x 10^"* cm-' 
(with r - 0.5 pc and Ar = 0.04 pc; see Figure [1]). Thus, they 
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Fig. 1. Emission structures from our 2D hydro simulation, for the same moments in time as panels 2 to 4 of Fi dA.4i i.e., 9000 yr, 
1 5 000 yr, and 1 8 000 yr after the onset of the BSG wind. Upper panel: Emissivity of the gas with 1 0"* < T < 1 0^ K, in the simulation 
plane. Lower panel: Projections of the 3D emission obtained by assuming rotational symmetry of the 2D structures of the upper 
panel, viewed with an inclination angle of 60°, constructed with a resolution of 400x400 pixels. 



shine with a total luminosity of 4 x 10^"^ erg s~' , i.e. roughly one 
solar luminosity, with a time scale of Ti-ad = l/u ^ 9000 yr, 
where u - jnkT is the internal energy of the gas, and T ^ 



10-23 gcm-3; 



10^ K and n ^ 13 cm"^ (corresponding to p 
Fig. IA.4b . The total radiated energy of the polar caps is about 
Emd - Ty^dL - 10"*^ erg. This corresponds well to the kinetic 
energy release due to the braking of the polar caps, which reach 
their maximum velocity of y ^ 130 kms"' at the time of colli- 
sion, where it is reduced to u ^ 50kms"' (Fig. |A.3T l. That is, 
AEkin = iAMAy2 ~ 8 x lO'^'^erg, with AM ^ 1.2 x IO^^m^ and 
Ay ^ 80kms '. This kinetic energy can be compared with the 
BSG wind kinetic energy, which, for M ^ 10"^'^ M© yr"' and 
f'wind - 300 km s"' (Fig. |A.2l l. yields ~ 1 .2 x 10"^^ erg over a time 
period of 9000 yr. Thus, the polar caps shine because the hour 
glass shaped BSG shell collides with the spherical RSG shell. 



A similar consideration could be made for the inner ring, 
which is produced by the collision of the BSG wind with the 
equatorial disk ejected by the central star during the phase of 
near critical rotation. The disk properties depend on the wind 
properties of the star during this phase. However, in particular 
their latitude dependence, can not be expected to be reliably pre- 
dicted within the current assumptions. The total mass of the disk 
is determined by the mass loss of the star at critical rotation. 






Fig. 2. Schematic representation of the interaction of the hour 
glass shaped BSG shell (blue) with the RSG shell (orange). The 
collision regions are marked with red color; they form the bright- 
est parts of the nebula. 



4. Discussion 

Figure|2]illustrates a simplified picture of the formation of mul- 
tiple ring nebulae, according to our model. It contains two kine- 
matic components: a stationary, spherical RSG shell and an ex- 
panding hour glass structure. The strongly emitting parts of the 
structure are the collision surfaces, which are marked in red in 
the figure. We believe that both kinematic components are re- 
alized by nature in the circumstellar medium of massive stars. 
RSG shells are unambiguously predicted (Garcia-Segura |2007] l 
and while they are not yet observationally confirmed, there 
seems to be no way to avoid their formation. Expanding hour 
glass structures, on the other hand, are a well documented fea- 
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ture in circumstellar nebulae of low and high mass stars (see 
Nota et al.'1995'; Langer et al. 1 19991 ) and are thought to be con- 
fined by a circumstellar disk in the equatorial plane of the central 
star 

A number of predictions emerge from this simple model. 
First, the collision starts about lO"* yr after the onset of the blue 
loop. This timing is set by the expansion speed of the BSG shell 
and the radius of the RSG shell. Second, the life time of the neb- 
ula is determined by the duration of the collision phase, as the 
emission time scale is shorter than that. In our example, this is 
about 10'* yr, or about 1% of the core helium burning life time. 
This provides an upper limit to the expected number of triple 
ring nebulae. Third, the rotation rate of the central star during the 
collision is high for a BSG, since it just about recovers from crit- 
ical rotation. At the time of maximum brightness of the nebula, 
the equatorial rotation rate of our central star is about 80kms"' 
(Fig. lA.2l i. Fourth, as all the material in the nebula is ejected after 
the first dredge-up phase of the central star, the nebula material 
is nitrogen-rich, here enhanced by a factor of 6.5, and carbon 
and oxygen depleted by factors of 6.5 and 1.5, respectively. We 
note that the level of N-enrichment predicted by current stel- 
lar evolution models is quite uncertain (see Hunter et al. 2008 ), 
but a RSG phase is still expected to produce some nitrogen en- 
hancement. Due to the assumptions of efficient rotational mix- 
ing (Heger & Langer 2000), the star and nebula in our model 
are more enriched than expected from non-rotating stellar mod- 
els. Fifth, one ingredient of our simple model, namely the RSG 
shell, is expected for massive stars, but not so for low mass star 
which produce planetary nebulae. Therefore, while quite analo- 
gous expanding hour glass structures are observed for both cases 
(Langer 2000), multiple ring nebulae formed through the colli- 
sion process shown in Fig.|2]are expected around massive stars, 
but not as planetary nebulae. In this sense, the polar caps ob- 
served around the blue supergiant Sher 25 might be considered 
as the first indirect empirical confirmation of a RSG shell. 

Previous models of multiple ring nebulae were mostly 
constructed in the context of the triple-ring structure ob- 
served aro und S N 1987A (Buitows et al. TgySl Crotts & 
Heathcote 120001 ). While single star models often fail to 
explain important features (e.g., Martin & Arnett 119951 
Meyer 119971 Woosley et al. 119971 ), many invoke rather com- 
plex binary phenomena (e.g., Podsiadlowski et al. 119911 
Blondin & Lundqvist [19931 Llyod et al. 1995; Morris & 
Podsiadlowski 20()5]). Whereas we do not attempt to reproduce 
the circumstellar medium of SN 1987 A, a single star approach 
with suitable choices for the major parameters in our model (ini- 
tial mass, initial rotation rate, metallicity) appears promising and 
will be pursued in the near future. The current failure of single 
star models to produce suitable blue loops and blue supergiant 
pre-supernova models may have to do more with missing physics 
in stellar evolution models rather than supporting the evidence 
for a binai-y progenitor of SN 1987A (Woosley et al. 1 19971 ) 

Various multiple ring nebulae around blue supergiants have 
been observed in the last 20 years (Smith et al. 2007). While 
our generic numerical model was not designed to correspond to 
any of them, many of the general properties of these nebulae 
are well reproduced. Most striking is the agreement of the emis- 
sion geometries. While the nebula around the B 1.5 la supergiant 
Sher 25 shows two polar caps and one equatorial ring (Brandner 
et al. 119971 ) and the other objects discussed by Smith et al. 120071 
rather show narrow rings, including the "twin" of the SN 1987 A 
nebula around HD 168625 (Smith 2007 ) all these structure occur 
as an evolutionary sequence in our model. Expansion velocities 
of the inner ring (~ 18 kms"') and the outer collision products 



(~ 50kms '), the spatial scale of about 1 pc, and the kinematic 
nebula age agree rather well with empirical values. The rotation 
velocity of our stellar model fits well to the derived value of 
~ 70kms ' for Sher 25 (Hendry et al. |2008] ). Central star and 
nebula of our model are nitrogen enriched, as are most of the 
observed nebulae. 

We note that the emission in our model is caused by com- 
pressional heating, which may be in conflict with evidence for 
photoionization being the dominant process in some observed 
multiple ring nebulae (see Smith et al. 2007V And indeed, look- 
ing at the density distributions shown in Fig. IA.4I which might 
resemble emission geometries in the pure photo-ionization case, 
the situation appears more complex. In our simulation, the thick 
RSG shell (Ar :^ r ^ 1 pc) collapses in two parts (at r ^ 1 .2 
pc and r ^ 1.7) due to a cooling instability. However, this col- 
lapse is questionable since it requires a long timescale — our 
shell has an age of close to 10^ yr, while in many cases the shell 
will be much younger at the time of collision — , and since the 
employed cooling function is uncertain for temperatures below 
10 000 K. Without this collapse, its density would be only about 
2 X 10"^^ g cm"-' (or 0. 1 particles/cm-'), which may render it un- 
observable even if it were photoionized. However, even in the 
case of the collapsed RSG shell as in our simulation, the colli- 
sion leads to a clear density enhancement. In panel 2 of Fig. |A.4| 
we see that in our model, the density enhancement in the po- 
lar caps is about a factor 5. This is to be considered a lower 
limit, as higher resolution models might approach the theoreti- 
cally expected enhancement factor of about 100, which follows 
from the (well realized) isothermal shock approximation and a 
Mach number of about 10. The lower panels of Fig. |A. 41 show 
that in order to represent the rings of SN 1987 A, further refine- 
ments are required, which, as we think, might be achieved by 
altering the properties of the RSG shell. For this particular case, 
this may indeed be justified, as the life time of the final RSG 
stage of the progenitor of SN 1987 A might have been quite short 
(Wooslev 1988 . Langer 1991a ). 

Despite that our model does not fit any of the observed cases 
in detail, the approximate agreement with most general proper- 
ties of this class of objects encourages to produce tailored mod- 
els for individual nebulae as next step. Our results indicate that 
stars with multiple ring nebulae might just have left the RSG 
branch — as stellar evolution models argued for the case of 
SN 1987A (Woosley 1988, Langer 1991a) and furthermore, that 
binarity may not be required to obtain multiple ring emission 
geometries. 
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Appendix A: Additional information 
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Fig. A.3. Blue supergiant shell velocity, covering a timespan of 
18 000 yrs, from its formation (f = 0), until its collision with the 
red supergiant shell. The time between two squares corresponds 
to 950 yrs. The time zero point corresponds to ~ 9000 yrs before 
the time of the first snapshot in FiglT] 



Fig.A.l. Evolution of our rotating 12 M© stellar evolution 
model in the Hertzsprung-Russell diagram, from the zero age 
main sequence to the pre-supernova stage. The initial metallic- 
ity of the model is Z = 0.02, and the initial equatorial rotation 
velocity is 328 km s"' . 



Fig. A.5. Movie description Movie of the projections of the 3D 
emission obtained by assuming rotational symmetry of the 2D 
structures obtained from the hydro simulation, viewed with an 
inclination angle of 60° (cf. Fig.[T] lower panel). It covers a time 
span of 17000 yr, starting at 2400 yrs after the time of the first 
panel of Fig. IA.4I The time difference between two frames is 
about 475 yr. 




t[yrx10^] 



t[yrx10^] 







■ 
1 


■ 






v^Dtm/c] 


■m 100 


J 


a 




3 




E 

>So 




^^^^::n:-:^ 


N^ 


n 






N . 



2.08 



1 
0.8 

0.6 

C 
0.4 

0.2 





Fig. A.2. Time dependence of various quantities of the employed 
12 Mq model, for about the final one million before core col- 
lapse. The upper panel shows stellar mass loss rate and terminal 
wind velocity. The lower panel depicts stellar rotational velocity, 
and the ratio of stellar and critical rotation rate, Q. 



S. M. Chija et al.: Multiple ring nebulae around blue supergiants. Online Material p 3 
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Fig. A.4. Snapshots of the density (left figures) and temperature (right figures) distributions of the circumstellar material around our 
12 Mq model after its first red supergiant stage. The first panel corresponds to beginning of the blue supergiant stage. The spherical 
red supergiant shell is situated between 1 pc and 2.2 pc. The hourglass shaped blue supergiant shell occupies the inner 0.4-0.5 pc. 
The second panel depicts the situation 4500 yr later, when the polar parts of the BSG shell collide with the RSG shell. Panels 3 
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